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QUANTITATIVE ANALYSIS OF CLAUSIUS INEQUALITY 
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Abstract. In the context of driven diffusive systems, for thermodynamic 
transformations over a large but finite time window, we derive an expansion of 
the energy balance. In particular, we characterize the transformations which 
minimize the energy dissipation and describe the optimal correction to the 
quasi-static limit. Surprisingly, in the case of transformations between homo¬ 
geneous equilibrium states of an ideal gas, the optimal transformation is a 
sequence of inhomogeneous equilibrium states. 


1. Introduction 

As discussed in thermodynamic textbooks, in a transformation between equilib¬ 
rium states a system necessarily goes through deviations from equilibrium which 
are small if the transformation is quasi-static. As clearly stated in Callen [5], 

A quasi-static process is thus defined in terms of a dense succession 
of equilibrium states. It is to be stressed that a quasi-static process 
therefore is an idealized concept, quite distinct from a real physical 
process, for a real process always involves nonequilibrium interme¬ 
diate states having no representation in the thermodynamic config¬ 
uration space. Furthermore, a quasi-static process, in contrast to a 
real process, does not involve considerations of rates, velocities or 
time. The quasi-static process simply is an ordered succession of 
equilibrium states, whereas a real process is a temporal succession 
of equilibrium and nonequilibrium states. 

Our aim is to develop the analysis, started in Hi, of real transformations for driven 
diffusive systems, both in equilibrium and nonequilibrium. As emphasized in the 
previous quotation this analysis will necessarily involve dynamical considerations 
that are outside the scope of classical thermodynamics. 

The dynamic evolution of driven diffusive system is described by the continuity 
equation together with the constitutive equation that expresses the local current 
as a function of the density and the driving field. The interaction with boundary 
reservoirs specifies the appropriate boundary conditions. A (real) transformation is 
thus defined by a choice of time-dependent driving field and chemical potentials of 
the reservoirs. Within this scheme, a dynamical derivation of the Clausius inequal¬ 
ity W > A F for isothermal transformations has been obtained in [2j S3 ■ Here W 
is the work done in the transformation and A F is the variation of the free energy. 
Moreover, we have shown that the Clausius inequality becomes an equality, i.e., 
W = A F in the limit of very slow transformations, that is in the quasi-static limit. 
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Since real transformations last a finite time, the quasi-static limit cannot be 
achieved. A meaningful issue is thus to describe the corrections to the quasi-static 
limit for transformations over a large but finite time window r. To describe the 
evolution of the system it is convenient to rescale time by introducing the dimen¬ 
sionless variable s = t/r where t is the original time variable. We then expand 
the evolution and the energy balance in powers of 1 /r and compute the first order 
corrections. 

Consider now transformations through equilibrium states namely, for which the 
stationary current, corresponding to the given external drivings at time s, van¬ 
ishes. In absence of an external held, transformations through equilibrium states 
are those in which the chemical potentials, while varying in time, are the same 
on each point of the boundary. For transformations through equilibrium states we 
show that, up to order 1/r 2 , W = A F + (1/r) B where B is a positive functional 
of the transformation, that is of the time dependent external drivings. For real but 
slow transformations we can thus optimize the dissipated energy by minimizing the 
functional B. Not surprisingly, we prove that for transformations between equilib¬ 
rium states (namely such that the initial and final states are equilibrium states), 
the functional B is minimized by transformations through equilibrium states. 

In the case of an ideal gas, the minimizer of B can be computed explicitly. 
Somehow surprisingly, for transformation between homogeneous equilibrium states 
(characterized by the absence of external field), the optimal transformation is a 
sequence of inhomogeneous equilibrium states. In other words, it is profitable 
to switch on an external field. In the context of Langevin dynamics, finite time 
refinements to the second law of thermodynamics have been discussed in [Tj, see 
also El for the case of jump Markov processes. 

For transformations between nonequilibrium states, the Clausius inequality W > 
A F does not carry any significant information. In fact, the energy dissipated along 
such transformations will necessarily include the contribution needed to maintain 
the nonequilibrium stationary states, which is infinite in an unbounded time win¬ 
dow. It is however possible to formulate a meaningful version of the Clausius 
inequality for nonequilibrium states by introducing a renormalized work W len that 
is defined by subtracting from the total work W the energy needed to maintain the 
nonequilibrium state. Within the setting of the macroscopic fluctuation theory [4] , 
a definition of renormalized work has been proposed in Ei following the point of 
view in [12] further developed in El HE]- The analysis of real thermodynamic trans¬ 
formations carried out in this paper and outlined above includes transformations 
between nonequilibrium states provided the work is replaced by the renormalized 
work. 

We draw the reader’s attention to the very recent paper [TO] - This paper analy¬ 
ses, in the context of Markov chains with finitely many degrees of freedom, problems 
similar to the ones discussed here. In particular, the authors introduce an optimiza¬ 
tion problem for the finite time correction to Clausius inequality, with motivations 
similar to ours. 


2. Clausius inequality and its nonequilibrium counterparts 

In this section we review the dynamical approach to thermodynamic transfor¬ 
mations introduced in [3 0 and developed in g]. 
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Hydrodynamical description. We introduce the hydrodynamic description of 
out of equilibrium driven diffusive systems which are characterized by conservation 
laws. We restrict to the case of a single conservation law, e.g., the conservation of 
the mass. 

We denote by A C the bounded region occupied by the system, by dA the 
boundary of A, by a; the macroscopic space coordinates and by t the macroscopic 
time. The system is in contact with boundary reservoirs, characterized by their 
chemical potential X(t,x), and under the action of an external field E(t,x). 

At the macroscopic level the system is completely described by the local density 
p{t,x ) and the local density current j(t,x). Their evolution is given by the conti¬ 
nuity equation together with the constitutive equation which expresses the current 
as a function of the density. Namely, 


f dtpit) + V -j(t) = 0, 
I j{t) = J(t,p(t)), 


( 2 . 1 ) 


where we omit the explicit dependence on the space variable x G A. For driven 
diffusive systems the constitutive equation takes the form 


J(t,p) = -D(p)Vp + x(p)E(t), 


( 2 . 2 ) 


where the diffusion coefficient D(p) and the mobility x{p) are assumed to be d x d 
symmetric and positive definite matrices. This holds in the context of stochastic 
lattice gases M- Equation (12.21) relies on the diffusive approximation and on the 
linear response to the external field. The evolution of the density is thus given by 
the driven diffusive equation 

dtpit) + V ■ ( X (p)E(t)) = V ■ (D(p)Vp) . (2.3) 

The transport coefficients D and y satisfy the local Einstein relation 

D(p) = x{ P ) f"(p)t (2-4) 

where / is the equilibrium free energy per unit volume. 

Equations m (Ei have to be supplemented by the appropriate boundary 
condition on dA due to the interaction with the external reservoirs. If A (t,x), 
x £ dA , is the chemical potential of the external reservoirs, the boundary condition 
reads 

(2.5) 


f'(p(t,x)) = xedA. 


If the chemical potential and external field do not depend on time, we denote by 
p = P\,e the stationary solution of (I2.3D . (12.51) . 


V • J{p) = V • ( - D(p)Vp + X (p) E) = 0, 
f'{p{x)) = A(x), x £ dA. 


( 2 . 6 ) 


We will assume that this stationary solution is unique. The stationary density 
profile p is characterized by the vanishing of the divergence of the associated current, 
V • J(p) = 0. A special situation is when the current itself vanishes, J(p) = 0; if this 
is the case we say that the system is in an equilibrium state; this can be viewed as 
a macroscopic counterpart to detailed balance. Conversely given a density profile 
p there is not a unique pair (A ,E) such that p = p\ t E- Indeed A is uniquely 
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determined by the second equation in (12.611 while the external field can be chosen 
in the form 

E = x~ 1 (p)(G + D(p)Vp), (2.7) 

where G is an arbitrary divergence free vector field. We note that for equilib¬ 
rium states there is a one-to-one correspondence between the pair (A, E) and the 
stationary solution of (12.611 . that is defined by choosing G = 0 in m- 

Homogeneous equilibrium states correspond to the case in which the external 
field vanishes and the chemical potential is constant. The stationary solution is then 
constant and satisfies f'(p\, o) = A. Inhomogeneous equilibrium states correspond 
to the case in which the external field is gradient, E = —VC/, and it is possible to 
choose the arbitrary constant in the definition of U such that U{x) = —X(x), x £ 
dA. By the Einstein relation (12.411 . the stationary solution satisfies —f'(p\ ,e{x)) = 
U(x ) and the stationary current vanishes, J(p\,e) = 0. Examples of inhomogeneous 
equilibrium states in presence of an external field are provided by a still atmosphere 
in the gravitational field or by sedimentation in a centrifuge. 


In this framework it is possible to define a thermodynamic functional V, called 
the quasi-potential, generalizing the free energy for systems out of equilibrium. It 
can be characterized as the maximal positive solution, vanishing when p = p, of 
the Hamilton-Jacobi equation, 

f SV dV f dV 

/ dxA7-T- ■ x(p)V "7 -/ dx-V-J(p)= 0. (2.8) 

J A °P op J A dp 

We now define the symmetric current Js by 

Js(p) = -x(p)V^. (2.9) 

dp 

Since the stationary density p is a minimum for V, then {dV/dp)(p) = 0. The 
symmetric current thus vanishes at the stationary profile, 

Jsifi) = 0. (2.10) 

We rewrite the hydrodynamic current as 

J(p) = J s (p) + Ja(p), (2.11) 

which defines the antisymmetric current Ja- 

In view of these definitions, equation (12.811 becomes 

[ dxJ s (p)-x(p)- 1 J A (p) = 0. (2.12) 

J A 

In the case of an equilibrium state the quasi-potential V = V\^e{p) is the local 
functional 

V x ,e{p)= [ dx(f(p)-f{p)-f(p)(p-p)), 

J A 

where p = p\ t E is the solution of (12.611 . 


(2.13) 
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Clausius inequality. The second law of thermodynamics can be expressed as 
follows. Consider a system in an equilibrium state in thermal contact with an 
environment at a given temperature. The system then undergoes an isothermal 
transformation to a final state. The Clausius inequality states 


W > A F, 


(2.14) 


where W is the mechanical work done on the system and A F is the difference 
of the free energy between the final and the initial state. If equality holds the 
transformation is said to be reversible. It can be implemented by performing very 
slow variations so that the system goes through a sequence of equilibrium states. 

We review the dynamical derivation of the Clausius inequality in n®. Consider 
a system in a time dependent environment, that is, E and A depend on time. The 
work done by the environment on the system in the time interval [0,T] is 


fC[o ,T] — 


[ dt | [ dxj(t) ■ E(t) - [ da X(t) j{t) • n|, (2.15) 

JO J A JdA ’ 

where n is the outer normal to d A and da is the surface measure on <9A. The first 
term on the right hand side is the energy provided by the external field while the 
second is the energy provided by the reservoirs. 

Fix time dependent paths A (t) of the chemical potential and E(t) of the driving 
field. Given a density profile po, let pit), j{t), t > 0, be the solution of (12.111 (12.511 
with initial condition po- By using the Einstein relation (12.4|1 and the boundary 
condition f'(p{t)) = A(f), an application of the divergence theorem yields 

W[o,t] = F{ P iT)) - E(p(0)) + [ dt f dx j{t) ■ (2.16) 

Jo J A 

where F is the equilibrium free energy functional, 


F{p) = [ dxf{p(x)). 
J A 


(2.17) 


Equation (12.1611 is not simply a rewriting of (12.151) . as it depends on a physical 
principle, the local Einstein relationship. 

Since the second term on the right hand side of (12.16|) is positive, we deduce the 
Clausius inequality (12.141) with A F = F{p\) — F{po) for arbitrary density profiles 
Po = p{ 0), pi = piT). Note that this derivation holds both for equilibrium and 
nonequilibrium systems. 

2.1. Renormalized work. The idea to define a renormalized work is to subtract 
the energy needed to maintain the system out of equilibrium. For time independent 
drivings, by the orthogonal decomposition (12.111) and (12.101) . J(p) = J a (p) is the 
macroscopic current in the stationary state. In view of the general formula for 
the total work (12.16L the amount of energy per unit time needed to maintain the 
system in the stationary profile p is 


dx J A (p) • x(p) Jf^ip)- 


(2.18) 


Fix now T > 0, a density profile po, and space-time dependent chemical poten¬ 
tials A it) and external field E(t), t £ [0,T]. Let (p(t),j(t)) be the corresponding 
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solution of (12.11) (12.51) with initial condition po. We define the renormalized work 
done by the reservoirs and the external field in the time interval [0,T] as 

Wfffa = W[ 0 ,t] ~ J q dt j^dxJ A (t,p{t)) • x(p(t)) _1 JA(i,p(t)) (2.19) 

where Ja(2, p) is the antisymmetric current for the system with the time indepen¬ 
dent external driving obtained by freezing the time dependent chemical potential A 
and external field E at time t. Observe that the definition of the renormalized work 
involves the antisymmetric current Jxit) computed not at density profile p\(t),E{t) 
but at the solution p(t) of the time dependent hydrodynamic equation. 

The definition (12.191) is natural within the macroscopic fluctuation theory and 
leads to a Clausius inequality. Indeed, in view of (12.161) and the orthogonality in 
(12.121) between the symmetric and the antisymmetric part of the current, 

W[ot] = F (P( T )) - F (Po) + J dt J^dx J S (t,p(t)) •x(p(t))- 1 J s (f,p(t)) ^ 2 20 ) 

> F(p(T)) - F(p 0 ). 

In the context of Langevin dynamics, a different definition of renormalized work 
has been proposed in (9], see [3] for a comparison. 

We obtain next a macroscopic version of the well known Hatano-Sasa inequality 
j6|. Consider the dissipation due to the symmetric current 

rT f f 5V 

J dt J dx J s {t,p{t))-x{p{t))~ 1 Js{t,p{t)) = -J dt J dx V-^- • J s (t,p{t)) > 0, 

where V* is the quasi-potential corresponding to the driving a = (A (t),E(t)) at 
frozen time t. Integrating by parts and using the orthogonality between J$ and J a 
we obtain 

J dt dx a • ^ > V t {p(T)) - V 0 (p( 0)) (2.21) 

If the initial state is the stationary profile for (A(0), 22(0)), then the right hand side 
is V t (p(T)) > 0. 


3. Finite time thermodynamics 


In this section we develop an approach to thermodynamic transformations which 
takes into account the fact that any real transformation lasts for a finite time. 
We consider transformations over an interval of time [0, t] and we discuss their 
asymptotic properties for large r. In particular, for slow transformation, we shall 
obtain the correction of order i to the equality W = AF which holds in the quasi¬ 
static limit. We finally discuss which transformations minimize such a correction. 


Slow transformations. To analyze transformations over the interval [0, r] it is 
convenient to introduce the dimensionless variable s = t/r. A protocol is defined 
by a choice of the external drivings E(s,x), x £ A, and X(s,x), x £ 9A, s £ [0,1]. 
The transformation is then realized by 


E T (t) 

A T W 


E (t/r ), 
A (t/r ), 


t £ [0,t], 


(3.1) 


The asymptotics in which we are interested is for r large compared to the typical 
relaxation time of the system, corresponding to slow transformations. Let p T {t) 
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and j T (t), 0 < t < t, be the solution to the hydrodynamic equations ()2.1I) . (12.211 
and m with the slow external field E T and chemical potential A r , that is 

r d tP T+v.j(t/T,pT(t)) = o, 

< j T (t) = J{t/r,p T {t)) (3.2) 

{ f'(P T (t))\ aA = * T (t) 

where we recall that J(t, p) = —D(p)'Vp + x(p) F (t). 

For s £ [0,1], let p(s ) be the unique stationary solution of the hydrodynamics 
with external field E(s) and chemical potential A(s). When r is large the solution 
(p T ,j T ) has an expansion of the type 

P T {rs) = p(s) + lr(s) + o(i), , . 

f(Ts) = J(s,p(s )) + iff(s) + o(i) . 

By (13.21) we get the corresponding linear evolution equations for the first order 
correction (r, g), 


d s p(s) + V ■ g(s) = 0 

g (s) = -V • (D(p(s))r(s)) + r(s)x'(p(s))E(s) (3.4) 

r(s, x) = 0, x £ dA 


where we use the notation that, for a matrix A = (ay(a;))” 3 - =1 , V • A is the vector 
with *-th coordinate d Xj aij(x). Note that the system (13.41) has the form of a 
Poisson equation for r(s). 

Recalling the definition (12.151) of the work, by evaluating the energy balance 
(12.161) along the transformation ( p T ,j T ), we obtain 


f {p t ( t )) ~ f {p t ( 0 )) = t [ ds [ dx j T (ts) ■ E(s) 

Jo J A 

— r / ds da\(s)j T (ts) ■ h ~t ds dx j T (ts) ■ x(p T ( rs )) 1 .7 r ( TS )- 
Jo Jd A Jo J A 


(3.5) 


In this equation p T ( 0) and p t (t) are the initial condition for the hydrodynamic 
equation (13.21) with the external drivings given in (13.11) and the corresponding value 
of the density at time r. Due to the finite relaxation time of the system, /9 r (0) and 
P t (t) are not the stationary density profiles associated to the drivings (A(0), E( 0)) 
and (A(l), E(l)). At order 1/r, the difference between p T ( 0) and p(0) is obtained 
by solving the equation (13.41) for s = 0 (here d s p(s) plays the role of a given source). 
The analogous statement holds for the difference between p T (r) and p( 1). Observe 
that in this formulation the value of the density at time 0 and t can be exchanged 
so that a slow transformation from the final to the initial state can be obtained by 
time reversal of the protocol. 

We can analyze the equation (13.51) at the different orders in 1/r, obtaining an 
identity for each order. Direct computations yield that at order r the right hand 
side of (13.51) vanishes. 

At order r° we get the first non trivial relationship, 


F(p(l)) — F(p(0)) = f ds f dx E(s) ■ g(s) — f ds f da\(s)g(s)-n 

Jo J A Jo Jd A 


+ [ ds [ dxr(s)J(s,p(s)) ■ (x 1 ) 1 (p(s))j(s, p(s)), 
Jo J A 


(3.6) 
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where we used that p T ( 0) = p(0) + 0(1 /t) and p T (r ) = p( 1) + 0(1 /t). The above 
relation connects the variation of the free energy to the first order corrections to the 
solutions of the hydrodynamic equations. We observe that, if we consider transfor¬ 
mations between two equilibrium states, the last term in the right hand side of (13.61) 
vanishes when the intermediate states are also of equilibrium so that J(s, p(s )) = 0. 
However the transformation can go through nonequilibrium intermediate states. 


Quantitative analysis of Clausius inequality. Consider the equation (12.201) 
which expresses the energy balance in the time interval [0,r]. Recall that the last 
term vanishes in the quasi-static limit. We now compute its asymptotics when r 
is large and for a slow transformation given, as in ED, in terms of a protocol 
(Hs),E(s)), s G [0,1]. 

Rewrite equation (12.201) for a slow transformation, 
W{^ t] -[F(p^(t))-F(p^(0))] 

= f dt f dx J s (t/T,p T (t)) ■ x(p T (t))~ 1 Js(t/r,p T (t)). 

JO J A 

Recalling (12.91) . the symmetric part of the current is 


T f \ ( \V7 SVx W,E(«)(P) [o Q\ 

Js(s,p) = -X(P)V - x - ( 3 - 8 

5p 

where Vx(s),e(s) is the quasi-potential associated to (A (s),E(s)) (we regard s here 
as a fixed parameter). In view of ( 13 . 31 ) . the symmetric current has the expansion 


Js(s,P t (ts)) = -^x(p(s))V (C s 1 r(s )) +0 (^). 
where Cj 1 is the linear operator with integral kernel 


Hence, 


C s 1 (x,y) 


6 2 V\ ( s ),e( s )(p(s)) 
Sp(x)Sp(y) 


TO ^ [ F(p T (r )) - F(p T ( 0))] = i B + 0(£) , 
where the excess functional B is: 

B = I ds [ dx\/(C~ 1 r(s)) ■x(p(s))^(C~ 1 r(s)). 
Jo J A 


(3.9) 


(3.10) 


(3.11) 


(3.12) 


For a transformation between and through equilibrium states, If,™] coincides with 
the total work W[ o, T ]. Hence, the inequality B > 0 is a restatement of the second 
principle of thermodynamics and (13.111) expresses a quantitative version of the 
Clausius inequality. Note that, in the limit r —>- oo, all protocols realize the equality 
W = A F. On the other hand, for finite time r, this identity cannot be achieved 
and we can select an optimal protocol by minimizing B. 


4. Optimal transformations between equilibrium states 

We consider, for simplicity, a system in one space dimension, in the domain A = 
[—1,1], with diffusion coefficient D(p) and mobility x(p)- Since p(s) is determined 
by (A (s),E(s)) through (12.fill , the excess (13.121) is a functional B = B(X,E) of the 
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protocol (A(s), E(s)), s £ [0,1]. In (13.121) r(s) = r\( s ^E(s)( x ) is obtained by solving 
the following Poisson equation, derived from (13.411 : 

d s p{s) = A (D(p(s))r(s)) - V ( X ' (p{s))E{s)r(s)) , . 

r(s, ±1) = 0 . { > 

Given an initial state (Ao,I?o) and a final state (Ai,i2i), we want to minimize 
the excess B(X,E) in (13.121) as a functional of the protocol, with the constraints 
(A(0),£(0)) = (A 0 , Eq) and (A(l),£(l)) = (Ai,£7i). 

This problem is already relevant when the initial and final states (Xq,Eq) and 
(Ai, £d) are equilibrium states. It appears reasonable that, in this case, an optimal 
protocol will pass through equilibrium states (A(s), E(s)) at every time s. Wc will 
show that this is indeed the case. Moreover, an optimal protocol can be obtained 
as follows. Solve the system of partial differential equations 

d s p(s) + |V(x(p(s))Vtt(s)) = 0 

d s 7r(s) + 3 X , (p(s))(Vtt(s)) 2 = 0 (4.2) 

0) = P\o,E 0 , P( 1) = P\i,Ei , 7r(s, ±1) = 0 , 
in the unknown p(s) = p(s,x) 1 Tt(s) = n(s,x), (s,x) £ [0,1] x A. Set 

E( S ) = ^l^-Vp( S ) , X(s, ±1) = f'{p(s, ±1)), (4.3) 

which corresponds to the choice G = 0 in (12.71) . Equation (14.31) defines a transfor¬ 
mation between and through equilibrium states, the corresponding minimal value 
of the excess functional is then given by 

B op t = j j ds J^dxx{p{s))(\7n(s)) 2 . (4.4) 

We emphasize that both the stationary equations (14.21) and the corresponding 
minimal excess B 0pt in (14.41) do not depend on the diffusion coefficient D. In fact, 
in the rescaled time s = and in the asymptotics r —> oo, the system relaxes 
instantaneously, and therefore the value of D becomes irrelevant. 

Remark on boundary conditions. Note that in the minimization of B we have not 
fixed the value of r(s) at s = 0 or 1, which corresponds, in terms of the unsealed 
time variable t, to fix the values of the initial and final density profiles p T (0) and 
p T (t) only at the order 1, and not at the order 4 . On the other hand, we claim 
that optimizing B with the added constraints r(0) = r(l) = 0 the infimum does 
not change. Indeed, we can consider a sequence of protocols that are constant in 
time in the time intervals [0,e] and [1 — er, 1], and close to the optimal protocol for 
s £ [e, 1 — e]. In the limit e —> 0 the corresponding value of B approaches 73 op t- As 
a consequence, we deduce that, for all protocols (A (s),E(s)), 

W[ 0 en T] - (E(pi) - F(po)) > ±B opt + 0(£), (4.5) 

provided that p T (0) = po and p T (t) = p\. The previous equation differs from (13.111) 
since there the variation of the free energy is F (p T (t)) — F (p T ( 0)). Equality in (14.51) 
can be achieved by the limiting procedure described above. 
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Change of variables. It will be convenient to perform a change of variables in 
the space of states. Given a state (A, E), we associate to it the pair density-current 
(p, J), where p = p\ t E is the stationary density profile defined by equation (USD , and 
J = —D(p)Vp+x{p)E is the corresponding stationary current. The correspondence 
(A, E) i-a (p, J) is one-to-one and the inverse map (p, J) e-t (A, E) is given by 

A(±l) = /'(p(±l)), E=^(D(p)Vp + J). (4.6) 

Observe that, since we are in one space dimension, J is constant in x. Under this 
change of variables, equilibrium states (A, E) correspond to elements (p, 0) with 
vanishing current. 

In the new variables, the quasi-potential V = U(p, J;p) becomes a functional 
on the set of density profiles p : A —> R + , depending parametrically on (p, J). For 
J = 0 it is the local functional (12.131) . 


V (p, 0; p) = f dx(f(p)~ f(p)-f'(p){p-p)) . (4.7) 

J A 

While for arbitrary current J £ R, the quasi-potential solves the Hamilton-Jacobi 
equation (12.81) . that in the present variables reads 


dx 


, SV(Pi J\p) 

5p 


v(|-(U(p, J;p)-U(p,0;p))) +-^-) =0, (4.8) 


J_ 

x(p )' 


where we used the Einstein relation (E3D- 

In the present variables, the excess functional (13.121) becomes 


B = 


3 2 V(p{s), J(s);p(s)) 




r(s,y)) 


where r(s) = r(p(s),d s p(s), J(s)',x) solves 

f d s p(s) = V (x{p(s)) V r ( s )) - X x (p(I)) r ( s )^i s 

\ r(s, ±1) =0. 


(4.9) 


(4.10) 


If the initial and final states are in equilibrium, then an optimal protocol consists 
of a family of equilibrium states (p(s), 0), s £ [0,1]. This will be shown by proving 
that the excess functional B in (14.91) satisfies 


SB(p, J) 
SJ(s) 


_ = 0 , 

j=o 


*€[ 0 , 1 ]. 


(4.11) 


Indeed, this condition guarantees that stationary paths (p(s), J(s)), s £ [0,1], of 
the excess functional B(p,J) can be obtained as (p(s),0), s £ [0,1], where p(s) 
is a stationary path for the functional B(p, 0). The proof of (14.111) is detailed in 
Appendix lAl 


Hamiltonian structure. For transformations between equilibrium states, in view 
of (14.111) . we can restrict the functional B to transformations through equilibrium 
states (p(s),0), s £ [0,1]. Under this assumption, the excess functional B (14.91) can 
be rewritten as 

B = Jo ^ J\ ^ ( V ( \(p(s )) r °(P( s )’P( s ))} ) • 


(4.12) 
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By introducing 

tt(s, x) = 7 r(p(s), p{s); x ) = -2 D ^ S,X ^ r 0 (p(s), p(s); x ) , (4.13) 

X(p(s,x)) 

the excess functional B can be written as 

B = \J q ds J A dx x(p(s)){Vn(s)) 2 . (4.14) 

and equation (14.101) translates to the following equation for 7r: 

dsP( s ) + ^ v (x(p(s))Vtt(s)^ = 0 , tt(s,±1) = 0. (4.15) 

In the form (14.141) . the excess functional B can be interpreted as the action 
functional associated to the Lagrangian 

£(P,P) = ^ j^ d xx{p{ x ))(yxK(p,p\x)) 2 . (4.16) 

The corresponding Hamiltonian is 

n{p,Tt) = sup | J^dx TTp — C(p, ^)| = ^ J^dxx{p)(y x 7f) 2 - (4-17) 

A straightforward computation shows that (14.21) are the Hamiltonian equations for 
(14.171) . Note that, apart for a factor (14.171) coincides with the Hamiltonian of 
the macroscopic fluctuation theory £U Sec.IVB] in the degenerate case D = 0 and 
E = 0. 


5. Explicit minimizers 


Optimal transformations through homogeneous equilibria. We start by 
discussing how the excess functional B can be minimized if we restrict it to trans¬ 
formations through homogeneous equilibrium states. Namely, we consider B in 
(I4.14I) - (I4.15I) as a functional on paths p(s), s G [0,1], constant in x (which corre¬ 
sponds to having zero external field: E(s) = 0). 

Within this setting equation (14.151) for 7r(s,x) = n(p(s), p(s); x) becomes 

Att(s, x) = -2 ^ 7r ( S) ±1) = 0 , (5.1) 

X(p(s)) 


whose solution is 

tt(s, x) = fy,, (1 - x 2 ). (5.2) 

X{p(s)) 

In view of (15.21) . the functional B (14.141) . restricted to homogeneous density proto¬ 
cols p(s) s G [0,1], becomes 


Letting 



(5.3) 



ds (5 s $(p(s))) 


2 


we have 


(5.4) 
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Hence, the minimizer of this functional is obtained when 

d s $(p{s)) = ^ S J = $(pi) - <P(p 0 ) ■ (5.5) 

VX'(p(s)) 

Thus the minimal excess (minimizing among the homogeneous protocols) is: 


-Dopt — ^ 


3 L 


<h(pi) - $(p 0 ) 


(5.6) 


The protocol (15.51) corresponds to the one obtained in d Eq.(18)] in the context 
of Markov processes with finitely many degrees of freedom. However, the spatial 
structure of our setting allows to find better protocols. In other words, the protocol 
(15.51) is not a minimizer of the excess functional (Id.141) without the constraint of 
transformations through homogeneous equilibrium states. Indeed, the function 
n(t,x) in (15.21) does not solve Hamiltonian equation (14.21) . In fact we get 


d s ir(s, x ) + ^X , (p(s))(V7t(s, x )) 2 = (p(s)) ( x (jf ( s )) ) ^ ^ ’ 

which does not vanish unless \ is constant (as in the so-called Ginzburg-Landau 
model). This means, in particular, that the optimal protocol will not be a sequence 
of homogeneous equilibrium states. In the case of ideal gases the actual minimizer 
will be found next. 


Ideal gas. In the case \{p) = Pi e -g-, f° r ideal gases, the Hamilton equation (14.21) 
reads 

f <9 s p(s) + iV(p(s)V7f(s))=0 

| d s n(s) + |(V7f(s)) 2 = 0 (5.7) 

[ P(0) = P\ 0 ,E 0 , M 1 ) = PXi,Er , 7f(s, ±1) = 0 . 

In particular, the second equation is decoupled and it admits solutions with sepa¬ 
rated variables. In the case p(0) = 0 and p(l) = p±, as can be checked by direct 
computations, the solution is as follows: 

7f(s, x) = ^(1 - |x|) 2 , p(s,x) =-6(\x\ + s - l)p!. (5.8) 

s s 

The corresponding minimal value of the excess functional is 

B 0 pt = gPi • (5-9) 

This should be compared to the minimal value of B through homogeneous equilibria 
(ESI), which in this case is |pi, giving a flat reduction of 75%. 

The interpretation of the solution (15.81) is the following. At time s = 0 + in¬ 
ject the required total mass 2pi at the endpoints of the domain, giving a positive 
contribution to the functional B. Then switch on the field E = D(p)^- : which is 
concentrated at the points x = ±(1 — s), so that the density profile p remains a 
step function at all times. Observe that the field E is opposite to the current, so 
the work done by the field is negative, and thus it gives a negative contribution to 
the excess functional B. 
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Ginzburg-Landau. This model has a constant mobility x{p) = c - I n this case 
equations equations m are linear and the solution is immediate 


p(x, s) = p\ 0 ,E 0 + s Oai.Ei - P\ 0 ,E 0 ) 
■k{s,x) = {p\ ltEl - P\o,Eo) (x 2 - 1) . 


(5.10) 


In particular, the optimal protocol is a sequence of homogeneous equilibrium states. 


6. Conclusions 


We have reviewed, in the context of driven diffusive systems, the macroscopic 
fluctuation theory approach to non-equilibrium stationary states [4]. In particular, 
we discussed the notion of renormalized work for which it is possible to prove 
a meaningful version of the Clausius inequality for transformation between non¬ 
equilibrium states. In the quasi-static limit this inequality becomes an equality. 

The main purpose of the present paper has been a quantitative discussion of the 
energy balance for real transformations, that is transformations lasting a large but 
finite time window. By rescaling the time variable we have obtained a relationship 
between the variation of the equilibrium free energy evaluated at the stationary 
density profiles and the first order correction to the hydrodynamic equation, see 
Eq. (1 3 6 [) . We have then introduced the excess functional B (see Eq. (13.1111 - (13.121) 1 
which accounts for the excess of the renormalized work with respect to the variation 
of free energy. Finally, we have analyzed the minimization of B. It is remarkable 
that, for transformations between homogeneous equilibrium states the optimal pro¬ 
tocol for B is not a sequence of homogeneous equilibrium states. This is due to the 
fact that in the framework of driven diffusive systems the spacial structure plays 
a nontrivial role and the system has a finite relaxation time. This result can be 
compared with the optimal protocol derived in El in the context of Markov pro¬ 
cesses with finitely many degrees of freedom. The optimal protocol for B has been 
explicitly computed in the case of an ideal gas in 1 space dimension, and it exhibit 
peculiar features. We also mention that, from a mathematical point of view, the 
optimization problem of B can be recast as a suitable optimal mass transportation 
problem. 

For real transformations between equilibrium states, the optimality criterion 
based on the minimization of the functional B in (13.121) appears a natural choice. 
On the other hand, for transformations between stationary non-equilibrium states, 
B only accounts for the excess of the renormalized work with respect to the variation 
of free energy. Therefore, a selection criterion based on the minimization of B in 
this case is meaningful when the subtracted counter-term f Q dt f A dx JA(t, p(t)) ■ 
X(p(t))~ 1 JA(t, p(t)) in (12.191) can be disregarded. This makes sense when the energy 
needed to maintain the stationary states is supplied by free unlimited sources, e.g. 
the solar energy. 


Appendix A. Analysis of the excess functional 


Preliminary computations. When the quasi-potential is local, i.e. it is as in 
(1X71) . then the operator C 1 in (13.101) is diagonal, and, in particular, it is given by 
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The quasi-potential V = V(p, J; p) has a minimun (equal to 0) in p = p, hence 
is has the general form: 

V(pJ',P) = \ [ dx f dyC~ 1 {p,J-,x,y)(p(x)-p{x)){p(y)-p(y))+o{p-p) 2 , (A.2) 
z J A J A 

where C(p , J]X,y) is the limiting covariance of density correlations. By setting 

the kernel C _1 in (IA.2I) satisfies 


(A.3) 


C 1 (p,J;x,y) = ^^-5(x-y) + JT(p-,x,y) + o(J). (A.4) 

X(p[x)) 

In order to prove (14.111) we need to find an equation for T(p; x, y). From (IA.2I) and 
(1A.4|) . we get 


6V{ f\ J } p) = 2j^(p(x) - P(x)) + J [ dy T(p; x, y)(p{y) - p(y)) + ■■■■ (A.5) 
Op{x) x{p{x)) J A 

We use equation (IA.5I) to expand the Hamilton-Jacobi equation (14.81) at order 2 in 
(p — p) and at order 1 in J. We get, after some simple algebraic manipulation and 
an integration by parts: 


J dx J^dy (p(x) - p(x))(p(y) - p(y))^^jjV x (. x(p(x))V x r(p ; x, y)J 
= jteWx) - - pM>) ■ 


(A.6) 


Since p(x) — p(x) is arbitrary, we conclude that T(p; x, y) (which is symmetric with 
respect to the exchange of x and y) satisfies 


(L x (p) + L y (p))T(p;x,y) = ~(R x (p) + R y {p))8(x - y ), 
r (p;aqy)|a(AxA) =0, 


where L(p) and R(p) are the differential operators defined by 

L m = ^(x(pw), W = f®v(A|y). (A.8) 

Equation (1A.7D is the desired equation for T(p;x,y). 


Recall that, given s € [0,1], the function r(s) : A —>• K defined by (14.10(1 . 
depending on the variable £ £ A, is a functional of p(s) : A —>• K + , d s p(s) : A —>• R, 
and J(s) £ R. Namely, we should denote r = r(p, p, J; x), and it is defined by the 
equation: 


v (* m (^l r ))- jv (m r ) =0 ’ r(±1)= ^ 

In the following we shall denote 


(A.9) 


r 0 (p,p;x) := r(p,p,0;x) 
We also introduce the new function 

dr{p,p, ■/; x) 


7(p,p;x) := 


dJ 


j=o 


In order to prove (14.111) we need to find an equation for q(p, p; x). 


(A.10) 

(A.11) 
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If we take the derivative of both sides of (IA.9I) with respect to J and we let 
J = 0, we get, after multiplying both sides by 

L (P ) (^jr^7 (p, P)) = R{p)ro(P, P ). 7(±1) = 0 , (A. 12) 

where L(p) and R(p) are the differential operators (|A. 8 |) . Equation (IA. 121) is the 
desired equation for 7 . 


Proof of (14.111) . By the definition (14.91) of the excess functional B and equations 
(1A.1I) . (IA.3I) . (1A. 101) and (IA.11H . we have 

SB(p, J ) 


SJ(s ) 


= 2 


J—0 


dx x(p(s, x)) (v x ( D ^ S ’ X ^ r 0 (p(s), p(s); x)) ) 
V Vx(p(s, x)) )) 


x (/ a dyV x T(p{s); x, y)r 0 (p{s), p(s);y) + V x l(p(s), p(s); x)) ) . 

XPS,X (A.13) 

Integrating by parts and recalling the definition (I A. 81) of the differential operator 
L(p), we can rewrite the right hand side of (IA.13I) . divided by —2, as 


j dx J dyr 0 (p(s) 1 p(s); x)r 0 (p(s), p(s); y)L x (p)r(p; x, y) 


+ / dxr 0 (p{s), p(s);x)L a 


X(p(s,x)) 


Using equation (1A.7[) . we get 


dx dyr 0 (p(s),p(s)-,x)r 0 {p(s),p(s)- : y)L x (p)r(p;x,y) 
J A J A 

= / dxr 0 (p{s),p{s)-,x)R x {p)r 0 (p{s),p(s)-,x). 

J A 

Moreover, using equation (IA.12I) . we get 


dx r 0 (p(s), p(s); x)L x (p) (- 


■X(p(s,x)) 

= / dxr 0 (p(s),p(s)-,x)R x (p)r 0 (p(s),p(s)-,x). 


(A.14) 


(A.15) 


(A.16) 


Combining (]A.14[) . (IA. 151) and (IA.16I) . we get (14.111) . 
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